The Falicov— Kimball model in external magnetic field: orbital effects 
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We study thermodynamic properties of the two-dimensional (2D) Falicov-Kimball model in the 
presence of external magnetic field perpendicular to the lattice. The field is taken into account 
by the Peierls substitution in the hopping term. In the non-interacting case the field dependent 
energy spectrum forms the famous Hofstadter butterfly. Our results indicate that for arbitrary 
nonzero interaction strength and arbitrary magnetic field there is a gap in the energy spectrum 
at sufficiently low temperature. The gap vanishes with increase of temperature for weak coupling, 
however, it persists at high temperatures if the coupling is strong enough. Numerical results have 
been obtained with the help of Monte Carlo technique based on a modified Metropolis algorithm. 



I. INTRODUCTION 

Strongly correlated electron systems have generated in- 
terest over the last few decades. It is widely accepted that 
many phenomena in condensed matter physics are con- 
nected with electronic correlations. Moreover an exper- 
imental progress made both in confining quantum gases 
and in preparation of nanosystems lead to requirement 
of theoretical models, that allow description of correlated 
quantum systems in periodic potentials and under an 
influence of the external magnetic field. Unfortunately, 
taking into account of all these three factors simultane- 
ously is intractable. While the problem of electrons in 
a periodic potential under influence of the external mag- 
netic field has been investigated since the beginning of 
the quantum mechanic, its solutions are known only in a 
few cases. In particular, one may prove that in two di- 
mensional fermionic gases orbital effects due to the mag- 
netic field lead the energy spectrum to form the famous 
Hofstadter butterfly. Incorporation of the electron corre- 
lations in such a model encounters yet unresolved prob- 
lems, mostly because of lack of relevant mathematical 
methods. 

The Hubbard model is frequently used as a starting 
point for studying strongly correlated systems. Origi- 
nally used for a description the metal-insulator transi- 
tion, it reveals interesting solutions, describing various 
phenomena observed in strongly correlated systems. De- 
spite the simplicity of the model, only a few rigorous 
results are known, mostly in one- or infinite-dimensional 
cases. Other results have been obtained with the help of 
approximate methods. Some attempts have been made 
to analyze the orbital effects in the Hubbard model in an 
external magnetic field. Some of them, based on an exact 
diagonalization, [l| were obtained for relatively small sys- 
tems and suffer from the finite size effects. The other at- 
tempts based on the mean-field approximation (MFA) [2| 
are questionable due to the limited applicability of the 
MFA in low dimensional systems. 

The lack of exact solutions for the Hubbard model 



and ambiguity of solutions obtained within approximate 
methods encouraged us to study even simpler model, i.e., 
the Falicov-Kimball model. It was proposed by Hubbard 
and Gutzwiller (for review, see, e.g., Q, as a simpli- 
fication of the Hubbard model and further was redevel- 
oped by Falicov and Kimball to study phase transitions 
in rare earths and transition metals. [5| The model is a 
limiting case of the asymmetric one-band Hubbard model 
where the mass of spin-down electrons goes to infinity. 
It describes a system consisting two kinds of fermions. 
One of them are itinerant particles and the others are 
massive and therefore localized. The only interaction in 
the Falicov-Kimball model is the on-site Coulomb re- 
pulsion between itinerant and localized particles. In the 
second quantization the model is described by the follow- 
ing Hamiltonian: 



H = 



(id) 



(1) 
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where t is the hopping integral, cj (/J) are the creation 
operators of an itinerant (localized) fermion at site i and 
U is the Coulomb interaction. 

While the Falicov-Kimball model is much simpler than 
its predecessor, it still cannot be rigorously solved in a 
general case. Fortunately, a significantly larger number of 
exact results is known for the the Falicov-Kimball model 
than for the Hubbard model. Q One of the most impor- 
tant theorem, proved by Kennedy and Lieb, [7j states 
that at low enough temperature there is a long range 
order for lattices of dimensionality greater than one for 
certain fillings and for all values of U. It has been shown 
that the ordering comes from effective correlations be- 
tween states of one kind, even though the model does not 
contain a direct interaction of this type. There are also 
important approximate results for the Falicov-Kimball 
model (see, e.g., Ref. §)• 

It has been demonstrated that in low dimensions the 
Falicov-Kimball model can effectively be analyzed with 
the help of the classical Monte Carlo (MC) method. 
[t| [T3| The method is general and it is possible to ap- 
ply it for different lattice geometries and for arbitrary 
fillings. The most important aspect is that the model 
allows for investigation of much larger systems than the 
original Hubbard model does. This is of crucial impor- 
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tance in the presence of magnetic field that results in a 
formation of the cyclotron orbits. It is clear that diame- 
ters of the orbits should not exceed the linear system size. 
Consequently, in the case of the Hubbard model one can 
investigate only extremely high magnetic fields when the 
magnetic flux through the lattice cell is of the order of 
the flux quantum. This limitation is significantly relaxed 
in the case of the Falicov-Kimball model. 

In the present paper we investigate a two-dimensional 
Falicov-Kimball model in a presence of perpendicular 
magnetic field. Since we analyze the spinless Falicov- 
Kimball model, the Zeeman term is absent. The effect of 
the Zeeman splitting was analyzed in, e.g., Refs. [lfj]. F° r 
U = the model reduces itself to the Azbel-Hofstadter 
model, the solutions of which form the famous Hofstadter 
butterfly. The solutions of the U ^ case give important 
guidelines how the electronic correlations modify the fine 
structure of Hofstadter butterfly. 

The outline of the paper is as follows. Section 2 briefly 
describes a model and a variation of the Monte Carlo 
method which is used to study systems with both clas- 
sical and quantum degrees of freedom. In Section 3 we 
present results obtained for Falicov-Kimball model in an 
external, perpendicular magnetic field. Section 4 con- 
tains summary and conclusions. 

II. MODEL AND COMPUTATIONAL METHOD 

In our study we analyze the extended Falicov-Kimball 
model, described by the following Hamiltonian: 

W = -£ ty (A)c\ Cj + uj2 fU4<*> ( 2 ) 

(i,J> 3 

where fy (A) is the hopping integral depending on the 
magnetic field through the Peierls phase factor: 

Uj (iS) = t exp ^ J ' A ■ dr^j , (3) 

and A = B(—ay, (I — a)x, 0) is the vector potential with 
the parameter a G [0, 1] that allows one to distinguish 
between the symmetric gauge (a = 1/2) and the Lan- 
dau gauge (a = 0). In the numerical calculations we 
have used the symmetric gauge. The same Hamiltonian 
was used by Gruber et al. to analyze flux phases in the 
Falicov-Kimball model. [l2| 

In our simulations we use a modified Metropolis 
algorithm. [loj As our system contains both itinerant 
fermions and localized particles, we use the grand canon- 
ical partition function in the following form: 

Z = Y J ^ee- m{c) ~^\ (4) 
c 

where C describes configuration of the localized states, /3 
is the inverse temperature and N is the operator of total 



number of itinerant fermions. For a given configuration 
C the Hamiltonian H.(C) can be diagonalized numerically 
and summation over fermionic degrees of freedom gives: 

z =En{ i+e ^" (cWi1 }' (6) 

C n 

where E n (C) is n-th eigenenergy of TL(C). Introducing 
the free energy of the mobile particles: 

^(C) = -i£ln{l + e-^ c )-"]} (6) 

n 

the partition function can be written in a form analogous 
to that used for the Ising model: 

Z = ^e-^ C \ (7) 
c 

where the difference from the Metropolis algorithm is 
that we use the electronic free energy instead of the in- 
ternal energy. MC simulations allow us to estimate the 
partition function and thermodynamic functions such as 
the specific heat. Position of the peak in the specific 
heat and in the charge density wave (CDW) susceptibil- 
ity allows one to determine the temperature T c of the 
transition between the ordered and disordered phases. @ 
In the present paper we restrict our considerations to 
a special case of half-filling both for itinerant and local- 
ized states. We have performed simulations of on square 
lattices with sizes from 10 x 10 to 30 x 30. Here, we 
present results for a 20 x 20 lattice. Although, for such a 
system we are unable to analyze the very fine details of 
the energy spectrum, it is possible to obtain the general 
structure of the Hofstadter butterfly within a reasonable 
computation time. 

III. RESULTS 

In order to the investigate thermodynamic properties 
of the Falicov-Kimball model with orbital effects due to 
magnetic field we performed simulations for different val- 
ues of the interaction strength U/t and for the entire 
range of the magnetic flux penetrating the lattice. It is 
convenient to use a dimensionless quantity a — 4>/4>o, 
where 4>q is the flux quantum and a € (0,1). It was 
shown by Kennedy and Lieb 0] that at low temperature 
the localized particles in the half-filled two-dimensional 
Falicov-Kimball model form a checkerboard pattern. At 
the same time the itinerant particles show the CDW or- 
dering. To reveal the most interesting properties of the 
model we have performed simulations for temperatures 
below, near and above the critical temperature T c of the 
transition between ordered and disordered phase. The 
critical temperature was determined from the peak in the 
specific heat versus temperature plot. The main results 
are presented in Fig. [TJ In the case of weak interactions, 
when U <t, the general structure of the Hofstadter but- 
terfly is nicely reproduced. At temperatures below T c the 
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FIG. 1: Hofstadter butterflies for different U/t and k B T/t. 
Figures were obtained for a lattice of 20 x 20 sites in half 
filling for both itinerant and localized particles. Darker color 
corresponds to higher value of the density of states. 
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FIG. 2: Examples of configurations of the localized particles 
at different temperatures in the absence of magnetic field. 



structure is similar to the energy spectrum of the nonin- 
teracting fermions, but is split near the Fermi energy by 
the band gap. With an increase of the temperature the 
gap diminishes and vanishes completely at the T c . Above 
T c the structure of the field dependent density of states 
looks similar to the Hofstadter butterfly, but is smeared 
because of averaging over disordered states. The effects 
of temperature on the configuration of the localized par- 
ticles is illustrated in Fig. O The presented configura- 
tions are "snapshots" of the evolution of the system dur- 
ing a MC run. Higher interaction strengths impact the 
energy structure more strongly. While for the interme- 
diate interaction strengths with U » At at temperatures 
below T c the energy spectrum still reproduces the Hof- 
stadter butterfly split by the band gap, increasing the 
temperature no longer reproduces the structure appear- 
ing for free fermion gas. Again at the temperature of T c 
the gap vanishes, but for higher temperatures the energy 
spectrum seems to be weakly dependent on the magnetic 



field. For a strong Coulomb interaction with U > 8t the 
gap persists at arbitrary temperature and its width is not 
smaller than U — 8t. An important fact is that at zero 
temperature in the absence of the interaction the mag- 
netic field switches the system between an insulating (or 
half metallic in an infinite system) and a metallic: phase, 
what results from the structure of the noninteracting Hof- 
stadter butterfly. In the presence of the interaction this 
behavior is not always reproduced. At temperatures be- 
low T c for all analyzed values of the ratio UJ t the sample 
is always in the insulating state, same as for U > 8t. 
On the other hand, when the temperature is above T c 
for weak and intermediate interaction U < 2t the mag- 
netic field switches the system from metallic to insulating 
phase (Fig. EJ). 



IV. CONCLUSIONS AND OUTLOOK 

We have shown that the numerical analysis of the 
Falicov-Kimball model in external magnetic field allows 
one to investigate the impact of the electron correlations 
on the Hofstadter butterfly. We investigated thermody- 
namics of the model with the Monte Carlo algorithm for 
various values of Coulomb interaction, magnetic flux and 
temperature. Simulations were carried out for lattices of 
10 x 10 to 30 x 30 sites for half filling for both kinds of par- 
ticles. For all probed values of U at temperatures below 
the T c states are insulating independently of magnetic 
field. For high values of the interaction strength U > 8t 
samples seem to be insulating at arbitrary temperature. 
An interesting result is that for weak and intermediate in- 
teraction strength the presence and magnitude of the en- 
ergy gap depends on the magnetic field in irregular man- 
ner, as in the noninteracting case. We show that for all 
investigated values of Coulomb repulsion there is a signif- 
icant smearing of the fine fractal structure of the energy 
spectrum. It occurs due to a broadening of quasiparti- 
cle levels. Nevertheless, it seems that the main branches 
of the Hofstadter butterfly survive (though split by the 
energy gap) in the presence of Coulomb interactions. 

The same analysis can be repeated for the Falicov- 
Kimball model away from half filling. In this case at low 
temperature the localized particles form patterns which 
arc incommensurate with the underlaying lattice. 
what results in a rich structure of the density of states 
even without magnetic field. The fractal structure of the 
Hofstadter butterfly is a result of an interplay between 
two length scales: lattice constant and the Landau ra- 
dius. In the case of the Falicov-Kimball model away from 
half filling there is an additional length scale, namely the 
period of the pattern formed by the localized particles. 
Therefore, one may expect the structure of the density 
of states versus the flux to be even richer than the Hofs- 
tadter butterfly. 

Another extension of the present approach would be to 
take into account the spins of the itinerant and/or local- 
ized particles. It would lead to the Zeeman splitting of 
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FIG. 3: Density of states obtained for various values of the magnetic flux a and Coulomb repulsion U ft. All densities of states 
were obtained from simulation at temperature T sa 1.5 T c , for a lattice of 10 x 10 sites with fixed boundary conditions. 



the energy levels, what in turn would result in drastic and 
nontrivial changes of the density of states at the Fermi 
level. [141] Therefore, in such a case the simple descrip- 
tion of the conditions for the metal-insulator transition 
presented in the previous section would not be valid any 
more. 
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